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Abstract:  This  research  project  has  addressed  critical  issues  concerned  with  the  mathematical 
modeling  of  phase  transitions  in  solids.  Detailed  microstructural  analysis  has  enabled  us  to  clarify 
the  role  of  stress  waves  in  the  propagation  of  phase  boundaries.  In  particular,  the  potential  for  phase 
boundary  movement  to  damp  out  stress  waves  has  been  demonstrated  for  certain  classes  of  material 
response.  Nucleation  of  new  phase  domains  by  highly  energetic  processes  has  been  successfully 
modeled  and  new  analytical  procedures  have  been  developed  for  the  predictive  response  of  phase 
transforming  media  during  high  energy  impact  Conversion  of  mechanical  energy  to  thermal  ener¬ 
gy  has  been  studied  by  means  of  an  extended  theory  which  incorporates  temperature  effects.  The 
role  of  these  dynamical  events  on  the  response  of  processes  and  devices  at  the  engineering  level 
points  to  the  utility  of  a  higher  order  mathematical  description  capable  of  capturing  cumulative  mi¬ 
crostructural  effects.  To  this  end  we  have  also  developed  mathematical  protocols  capable  of  track¬ 
ing  the  evolution  of  thermomechanical  austenite/martensite  phase  variants  due  to  generalized  con¬ 
ditions  of  loading  and  heating.  The  associated  mathematical  model  is  capable  of  capturing  super¬ 
elastic  response  and  two-way  shape  memory  in  diermoelastic  materials.  Ohmic  heating/convecdve 
cooling  is  a  lUcely  activation  mechanism  for  smart  devices  utilizing  these  materials,  and  our  numer¬ 
ical  simulations  have  successfully  replicated  these  processes  including  a  smooth  transition  from 
isothermal  to  adiabatic  response  as  the  loading  rate  is  increased  in  a  heat  convective  environment 
The  effect  of  incorporating  these  highly  nonlinear  materials  in  composite  layered  arrangements  and 
thin  frlm  structures  has  also  been  examined.  This  research  project  has,  to  date,  generated  eleven 
archival  papers  and  conference  proceedings,  contributed  to  eighteen  technical  presentations,  and 
supported  two  Ph.D.  graduate  students. 
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1.  Introduction 

The  puipose  of  this  research  project,  as  stated  in  the  original  research  proposal,  was  to  address 
"critical  issues  concerned  with  the  mathematical  modeling  of  phase  transitions  in  solids.  The  spe¬ 
cific  phase  transitions  of  interest  are  those  which  are  rapidly  activated  by  either  mechanical  or 
thermal  means.  The  influence  of  stress  waves  will  then  exert  a  fundamental  role  on  the  triggering 
and  subsequent  propagation  of  phase  boundaries.  At  present,  the  ability  to  predict  the  dynamical 
response  of  phase  tran^orming  materials  is  insi^ficiently  precise  for  confidence  in  understanding 
how  these  materials  will  react  under  extreme  conditions  or  perform  in  novel  design  settings.  This 
research  is  designed  to  develop  analytical  methods  for  adtbressing  these  issues.”  The  original  re¬ 
search  proposal  then  proceeded  to  describe  in  more  detail  the  specific  problems  to  be  examined  and 
the  medio^  proposed  for  their  treatment 

We  are  happy  to  report  that  these  goals  have  been  met  Substantial  progress  has  been  achieved  on 
all  of  the  problem  areas  proposed  for  original  study.  This  is  evidenced  by  the  numerous  publica¬ 
tions  and  technical  presentations  which  have  resulted  from  this  research.  Two  graduate  students. 
Jack  Lin  and  Rtdph  Worthington,  have  benefitted  from  intense  involvement  in  this  research.  It  is 
anticipated  that  Jack  will  receive  Us  Ph.D.  in  May  1993  and  that  Ralph  will  receive  his  Ph.D.  some¬ 
time  in  1994.  In  addition,  our  accomplishments  have  suggested  new  avenues  of  inquiry,  some  of 
which  we  have  also  begun  to  pursue.  In  particular,  we  have  initiated  collaboration  with  an  industrial 
research  partner  which  has  made  possible  a  much  faster  degree  of  technical  t  ^nsfer  than  that  which 
would  be  possible  by  the  traditional  publication  &  presentation  process  (this  industrial  collabora¬ 
tion  has  also  resulted  to  date,  in  two  soon-to-appear  publications).  The  support  of  ARO  has  also 
permitted  us  to  obtain  a  certain  amount  of  leveraged  research  frinding  from  the  State  of  Michigan 
Research  Excellence  Fund  ($24,991,  10/89-9/90),  ($24,000,  10/90-9/91),  ($27,000,  10/91-9/92), 
($31,911, 10/92-6/93)  to  extend  some  of  the  issues  addressed  in  our  ARO  study  to  the  domain  of 
composite  material  structures. 

The  support  provided  by  this  research  grant  has  been  instrumental  to  the  success  of  the  principal 
investigator’s  academic  career  during  the  past  three  years.  During  this  time  1 1  papers  sup^rted  by 
ARO  have  either  appeared,  or  have  b^n  accepted  for  publication,  in  the  journals  of  reconl  for  me¬ 
chanics  and  applied  mathematics  (including  the  Archive  for  Rational  Mechanics  and  Analysis,  the 
Journal  of  Elasticity,  and  the  Intematioruil  Journal  of  Engineering  Science),  and  additional  papers 
are  in  various  stages  of  preparation.  Eighteen  presentations  have  ^n  delivered  at  external  confer¬ 
ences  and  other  universities.  The  principal  investigator  was  promoted  by  his  institution  from  As¬ 
sistant  Professor  to  Associate  Professor  with  tenure,  and  was  ^so  awarded  the  College  of  Engineer¬ 
ing  Withrow  Teaching  Excellence  Award  in  1991,  the  inaugural  year  for  this  award.  In  1992,  Ph.D. 
graduate  student  Jack  Lin  was  recognized  by  the  College  of  Engineering  as  the  outstanding  gradu¬ 
ate  student  in  the  Department  of  Materials  Science  and  Mechanics. 

During  the  period  of  time  encompassed  by  this  research  (8-1-89  to  12-15-92)  we  have,  of  course, 
continued  to  monitor  the  research  advances  obtained  by  other  investigators  in  this  area  (some  of 
whom  are  also  supported  by  ARO),  and  our  research  has  benefitted  from  this  scholarly  interaction. 
Now  that  the  success  of  our  mathematical  modeling  procedures  has  been  established,  we  are  be¬ 
ginning  to  seek  avenues  for  the  additional  application  of  our  results.  To  this  end  we  have  begun 
collaboration  with  materials  scientists  whose  specialities  are  in  materials  fabrication  and  process¬ 
ing.  In  particular,  as  of  this  writing,  research  proposals  are  pending  to  hopefully  enable  us  to  pursue 
such  extensions  (a  research  proposal  entitled  Mathematical  Modeling  arid  Characterization  of  Self- 
Biasing  Thermoelastic  Thin  Film  Microlaminates  on  High  Temperature  Polymeric  Substrates,  co¬ 
authored  with  materials  scientist  D.S.  Grummon,  has  bwn  submitted  to  the  ARO’s  URI-RIP  93 
program  (control  number  SM- 13548)). 

In  the  next  two  sections  of  this  final  report  we  list  the  publications  that  have  already  resulted  from 
this  research  project,  as  well  as  the  external  presentations  delivered  during  the  granting  period. 
Then  in  the  next  five  sections  we  summarize  in  more  detail  our  various  technical  achievements.  In 
this  description  we  underscore  not  only  how  we  have  fulfilled  the  goals  of  the  originally  proposed 
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project,  but  also  how  new  results,  both  by  us  and  others,  have  been  successfully  incorporated  into 
our  research  program.  We  also  touch  upon  some  of  the  new  avenues  of  opportunity  opened  by  our 
endeavors. 


2.  Publications  Acknowledging  ARO  Support 

The  following  eleven  papers  all  acknowledge  ARO  support  and  have  either  appeared  or  been  ac¬ 
cepted  for  publication.  These  papers  have  been  submitted  to  the  ARO  library  as  per  the  grant  spec¬ 
ifications  (except  for  the  Transactions  of  the  Army  Conferences  which  originated  from  ARO  pub¬ 
lication  sources),  hi  addition,  other  manuscripts  are  presendy  in  various  stages  of  preparation 
which  will  also  acknowledge  ARO  support 

Ivdiin,  Y.  and  T.  J.  Pence  (1992a).  A  constitutive  model  for  hysteietic  phase  transition  behavior,  (accepted 
for  publication  in  Int.  J.  Eng.  ScL). 

Ivshin,  Y.  and  T.  J.  Pence(1992b).  A  simple  mathematical  model  of  two-way  memory  effect  (to  appear  in 
Proceedings  of  die  1992  Int.  Cortf.  on  Martensitic  Transformations  (ICOMAT-^)). 

Lin,  J.  and  T.  J.  Pence  (1991).  Energy  dissipation  in  an  elastic  material  containing  a  mobile  phase  boundary 
subjected  to  concurrent  dynamic  pulses.  Transactions  of  the  Ninth  Aimv  Conference  on  Applied  Math¬ 
ematics  and  Computing.  F.  Dressel  ed.,  437-4S0. 

Lin,  J.  and  T.  J.  Pence  (1992a).  On  the  dissipation  due  to  wave  ringing  in  non-elliptic  elastic  materials,  (ac¬ 
cepted  for  publication  in  J,  Nonlinear  Sci.). 

Lin,  J.  and  T.  J.  Pence  (1992b),  Kinetically  driven  elastic  phase  boundary  motion  activated  by  concunem 
dynamic  pulses,  to  appear  in  Transactions  of  the  Tenth  Annv  Conference  on  Applied  Mathematics  and 
Computing.  F.  Dressel  ed. 

Pence,  T.J.  (1990a).  Phase  transitions  and  maximally  dissipative  dynamic  solutions  in  the  Riemann  problem 
for  impact  Transactions  of  the  Eighth  Army  Confer^  on  Anidied  Mathematics  and  Computing.  F. 
Dressel  ed.  817-832. 

Pence,  T.J.  (1990b).  An  inverse  Riemann  proUem  for  impact  involving  phase  transitions.  Transactions  of 
the  Eighth  Army  Conference  on  Applied  Mathematics  and  Computing.  F.  Dressel  ed.  833-846. 

Pence,  T.J.  (1991b).  On  the  erKounter  of  an  acoustic  shear  pulse  with  a  phase  boundary  in  an  elastic  mate¬ 
rial:  energy  and  dissipation.  J.  Elasticity  26, 95-146. 

Pence,  TJ.  (1992a).  On  the  mechanical  dissipation  of  solutions  to  the  Riemarm  problem  for  impact  involv¬ 
ing  a  two-phase  elastic  material.  Arch.  Ration.  Mech.  &  Anal.  117, 1-52. 

Pence,  TJ.  (1992b).  The  dissipation  topography  associated  with  solutitms  to  a  Riemann  problem  involving 
elastic  materials  undergoing  phase  transitions,  to  appear  in  Shock  Induced  Transitions  and  Phase 
Structures  in  General  Media,  R.  Fosdick,  E.  Durm  &  M.  Slemrod  eds.,  IMA  volume  series,  Springer- 
Verlag. 

Song,  J.  and  T.  J.  Pence  (1992),  On  the  design  of  three-ply  nonlinearly  elastic  composite  plates  with  optimal 
buckling  performance.  Structural  Optimization  5  (1992)  45-54 


3,  External  Presentations  During  the  Grant  Period 

Eighteen  presentations  were  given  at  technical  conferences,  or  delivered  during  invited  seminars 
at  other  universities: 


XXI  Midwestern  Mechanics  conference.  August  13-16, 1989.  Houghton,  MI. 
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SES  26"'  Annual  Meeting,  September  18-20, 1989,  Ann  Arbn,  MI. 

Invited  seminar,  Departmoit  of  Mechanical  Engineering,  The  Johns  Hc^dns  University,  Oct.  12, 1989. 

Invited  seminar.  Department  of  Applied  Mathematics,  University  of  Virginia-Charlottesville,  Oct.  13, 1989. 

Society  for  Natural  Philos(q;)hy  34"'  Meeting,  ^>ril  6-8, 1990,  Lincoln,  NB. 

Eleventh  National  Congress  of  ^rplied  Mechanics,  (with  J.  Song),  May  21-25, 1990,  Tucson,  AZ. 

Eleventh  National  Congress  of  ^rplied  Mechanics,  (widr  Y.  Ivshin),  May  21-25, 1990,  Tucson,  AZ. 

American  Society  of  Composites  5"*  Technical  Conference,  (with  J.  Song),  June  12-14, 1990,  East  Lansing,  MI. 
Eighth  Army  Conference  on  ^n>Ued  Mathematics  and  Computing,  June  19-22, 1990,  Ithaca,  NY. 

IMA  woricshop:  Shock  Induced  Transitions  and  Phase  Structures  in  Gerteral  Media,  Oct  15-19, 1990,  Miruieapolis. 
Invited  lecture,  Dqrt  of  Aerospace  Engineering  and  Applied  Mechanics,  University  of  Minnesota,  Oct  30. 1990. 
One  hundred  eleventh  Winter  Annual  Meeting  of  ASME,  Nov.  25-30, 1990.  Dallas.  TX. 

1991  Mathematica  Conference,  (with  Y.  Ivshin),  January  12-15, 1991,  San  Francisco,  CA. 

Invited  seminar,  Depts.  of  Civil  Engineering  and  Mechanical  Engineering,  Northwestern  University,  May  17, 1991. 
Ninth  Army  Conference  on  Applied  Mathematics  and  Computing,  June  18-21, 1991,  Minneapolis,  MN. 

SES  28*^  Armual  Meeting,  November  6-8, 1991,  Gainesville,  FL. 

Tenth  Army  Conference  on  Applied  Mathematics  and  (Computing  (with  J.  Lin),  June  16-19, 1992,  West  Point  NY. 
Int  Conference  on  Martensitic  Transformations,  (with  Y.  Ivshin),  July  20-24, 1992,  Monterey,  CA. 


4.  Phase  Boundary  Propagation  in  the  Presoice  of  Stress  Waves 

As  background,  phase  transformations  in  solids  can  be  modeled  in  a  purely  equilibrium  setting 
within  the  finite  deformation  theory  of  elasticity  provided  that  the  elastic  potential  energy  density 
of  the  material  admits  a  multiple  well  structure,  lliis  allows  the  governing  equations  to  lose  ellip- 
ticity,  consequently  the  solutions  of  boundary  value  problems  may  then  contain  internal  surfaces 
of  strain  discontinuity  that  separate  distinct  phase  regions.  Phase  boundaries  of  this  type,  and  the 
attendant  configurations  involving  co-existing  solid  phases  in  equilibrium,  have  been  studied  by 
numerous  authors,  including  Knowles  and  Sternberg  (1978),  James  (1979,  1981,  1983,  1986), 
Kikuchi  and  Triantafyllidis  (1982),  Fosdick  and  Macsithigh  (1983),  Gurtin  (1983),  Abeyaratne  and 
Knowles  (1987a,  b,  1988a,  b).  Silling  (1989)  and  Jiang  (1991).  A  consistent  feature  of  these  studies 
is  that  boundary  value  problems  which  would  be  well-posed  in  the  conventional  “elliptic  theory” 
(viz.  Ball  1977)  often  become  underposed  in  the  “non-elliptic  theory”:  namely  a  massive  loss  of 
uniqueness  prevails  if  only  the  usual  constitutive  and  mechanical  equations  of  equilibrium  are  re¬ 
quired  to  hold. 

These  different  equilibrium  configurations  will  typically  each  store  a  different  amount  of  mechan¬ 
ical  energy.  In  certain  cases,  configurations  that  are  global  energy'  minimizers  may  be  identified, 
although  it  also  may  be  the  case  that  more  than  one  configuration  provides  a  global  energy  mini¬ 
mum  state  (Ericksen  1975).  In  other  cases,  energy  minimal  equilibrium  configurations  might  not 
exist,  in  the  sense  that  a  lower  bound  energy  level  is  approachable  but  not  attainable.  This  may  ex¬ 
plain  the  phenomena  of  “frustration”  in  magnetic  materials  (James  and  Kindeilehrer  1990).  In  ad¬ 
dition,  when  the  governing  energy  densities  explicitly  account  for  the  crystallographic  symmetries 
of  the  competing  phases,  the  microstructural  “quest”  for  an  unattainable  minimum  energy  state 
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might  then  be  responsible  for  fme  variations  in  crystallographic  domain  structure  (Ball  and  James 
1987,  1990).  Energy  based  methods  for  determining  optimal  microstructural  types  and  arrange¬ 
ments  of  phases,  voids  and  other  defects  are  presently  an  active  area  of  study  (e.g.  Collins,  Kinder- 
lehrer  and  Luskin  1991;  Collins  and  Luskin  1991;  James  and  Spector  19^;  Kohn  1989,  1991; 
Bhattacharya  1991;  Fosdick  and  Zhang  1992).  The  presence  of  hysteresis  and  generalized  “mem¬ 
ory”  effects  in  the  quasi-static  evolution  of  equilibrium  states  suggests  that  energy  minimization 
might  not  be  the  sole  governing  agent  for  the  determination  of  equilibrium.  It  has  been  suggested 
(Abeyaratne  aiKi  Knowles  1988a)  that  kinetic  balance  relations  operating  at  phase  boundary  inter¬ 
faces  may  then  control  the  evolution  of  equilibrium  states. 

If  the  dynamics  of  these  systems  is  purely  elastic,  then  the  governing  field  equations  and  disconti¬ 
nuity  conditions  involve  only  the  addition  of  inertial  terms.  Phase  boundary  motion  is  possible  in 
this  setting  and  the  underposedness  problem  remains  (James  1980;  Shearer  1982,  1986,  1988; 
Slemrod  1983, 1989;  Hattori  1986a,b;  Pence  1986, 1987;  Traskinovsky  1987, 1990).  In  fact  if  one 
considers  a  body  in  equilibrium  which  is  then  dynamically  disturbed,  there  are  now  two  possible 
sources  of  nonuniqueness:  nonuniqueness  in  initial  conditions  stemming  firom  those  issues  that  are 
resident  in  the  equilibrium  theory  as  mentioned  above,  and  an  additional  purely  dynamical  nonu¬ 
niqueness  due  to  the  fact  that  the  dynamical  motion  from  a  fixed  set  of  initial  conditions  may  also 
admit  multiple  solutions  for  t  >  0. 

Turning  now  to  dynamical  nonuniqueness,  a  moving  phase  boundary  in  one  space  dimension  is  in 
general  either  a  point  source  or  a  point  sink  for  mechanical  energy.  Standard  interpretations  of  the 
second  law  of  thermodynamics  then  require  that  only  sink-like,  or  dissipative,  motion  is  admiss- 
able  (Dafermos  1973;  Abeyarame  and  Itoowles  1990a;  Gurtin  1991).  This  restriction  is  typically 
not  sufficient  to  wholly  resolve  issues  of  dynamical  uniqueness,  and  additional  criteria  have  been 
proposed  for  this  purpose  including:  requirements  of  maximizing  energy  dissipation  (Dafermos 
1973, 1989;  James  1980;  Hattori  1986a),  and  the  application  dynamical  kinetic  relations  governing 
phase  boundary  mobility  (Truskinovsky  1987;  Gurtin  and  Struthers  1990;  Abeyarame  and 
Knowles  1991). 

Dynamical  nonuniqueness  is  not  at  issue  in  the  event  that  the  equations  of  motion  involve  not  only 
the  inertial  terms  of  the  purely  elastic  theory,  but  also  an  additional  viscosity  term  (Slemrod  1983; 
Pego  1987;  Ball  et  al.  1991).  The  equilibrium  theory  is  not  affected  and  the  regularization  provid¬ 
ed  by  the  higher  order  derivative  delivers  a  unique  solution  to  problems  that  begin  from  a  fixed  set 
of  initial  conditions.  In  addition,  viscous  dissipation  takes  place  throughout  the  field.  One  might 
thus  expect  that  these  dynamical  solutions  would  then  have  two  dissipative  mechanisms;  viscous 
dissipation  throughout  the  field  and  localized  energy  sinks  provided  by  moving  phase  boundaries. 
In  fact  there  is  only  the  viscous  field  dissipation,  as  the  viscous  regularization  neither  permits  the 
motion  of  initial  phase  boundaries  nor  allows  for  creation  of  new  phase  boundaries  (Proposition  35 
and  Theorem  3.6  of  Pego  1987;  Theorem  4.10  of  Ball  et  al.  1991).  Modeling  of  other  physical  effects,  such 
as  capillarity  (Hagan  and  Slemrod  1983)  or  a  short  thermal  memory  (Niezgodka  and  Sprekels 
1988;  Niezgodka,  Zheng  and  Sprekels  1988)  can  give  rise  to  additional  higher  order  derivative 
terms  in  the  governing  equations;  the  solutions  of  these  equations  may  then  be  sensitive  to  various 
scalings  of  these  terms  even  when  their  individual  contributions  are  vanishingly  small  (Abeyarame 
and  Knowles  1990b). 

This  drastic  suppression  of  phase  boundary  motion  in  models  which  utilize  viscosity  terms  in  the 
governing  equations  of  motion  is  at  odds  with  the  observed  character  of  diffusionless  phase  trans- 
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fonnations.  The  austenite-martensite  phase  boundary,  for  example,  is  a  sharp  interface  which  is 
-  -  found  on  occasion  to  travel  at  speeds  approaching  the  shear  wave  speed  of  the  material  (Bunshak 
and  Mehl  1952;  Grujicic,  Olson  and  Owen  1985).  This  supports  the  consideration  of  dynamical 
models  for  phase  transitions  that  involve  elastic  dynamics  supplemented  by  appropriate  local  cri¬ 
teria  for  phase  boundary  mobility.  In  Pence  (1991a)  we  utilized  this  framework  to  investigate  the 
role  of  shear  waves  in  driving  a  phase  boundary  through  a  two-phase  material.  Specifically,  we 
considered  the  interaction  of  a  finite  amplitude  acoustic  shear  pulse  with  a  pre-existing  phase 
boundary  in  a  material  that  was  in  equilibrium  prior  to  the  introduction  of  the  pulse.  In  general, 
from  this  encounter  there  eventually  emerges  both  a  reflected  pulse  and  a  transmitted  pulse.  The 
dynamical  fields  associated  with  the  different  possibilities  for  phase  boundary  motion  were  com¬ 
pletely  parametrized  in  terms  of  the  phase  boundary  speed.  For  certain  ranges  of  iiutiating  shear 
pulse  amplitude,  a  completely  reflecting  solution  was  found,  while  for  an  in  general  different  range 
of  shear  pulse  amplitudes  a  completely  transmitting  solution  was  found,  llie  properties  of  these 
different  solutions  suggest  that  the  notion  of  phase  boundary  impedance  may  be  useful  in  deter- 
mirung  criteria  for  phase  boundary  mobility. 

Subsequently  (Pence  1991b)  these  results  were  examined  from  the  perspective  of  energy  and  dis¬ 
sipation.  It  was  shown  that  there  exists  at  most  two  motions  which  involve  no  dissipation  (corre¬ 
sponding  to  conservation  of  mechaiucal  energy),  that  there  exists  a  third  motion  which  maximizes 
the  mechanical  energy  dissipation  rate,  and  that  there  exists  yet  a  fourth  motion  which  maximizes 
the  total  energy  loss  during  the  encounter  (the  distinction  between  the  third  and  fourth  motions  is 
due  to  the  fact  that  certain  encounters  with  lower  dissipation  rate  have  longer  duration).  This  raises 
the  issue  of  whether  or  not  the  disturbance  applied  to  an  energetically  stable  equilibrium  state  pro¬ 
vided  by  such  a  finite  amplitude  shear  pulse  can,  through  subsequent  shear  wave  reverberations, 
completely  damp  itself  out,  thus  allowing  the  material  to  assume  a  new  energetically  stable  mini¬ 
mum  state  consistent  with  the  external  load  on  the  system  after  the  pulse  has  been  generated.  Non- 
dissipative  motion  does  not  allow  for  this  possibility.  However  it  was  established  (Pence  1991b) 
that  this  possibility  is  not  excluded  for  either  the  third  (maximally  dissipative)  motion,  or  the  fourth 
(maximal  loss)  motion.  In  fact,  the  optimal  energy  loss  for  a  single  maximally  dissipative  encoun¬ 
ter  is  70.5%  of  that  needed  to  attain  the  new  ground  state,  and  this  occurs  within  materials  for  which 
the  transformation  strain  is  96.7%  of  the  initiating  pulse  amplitude.  These  results  put  rigorous 
bounds  on  the  possibilities  for  energy  absorption  via  a  single  (and  hence  rapid)  phase  transforma¬ 
tion  event  The  energetic  behavior  of  a  phase  boundary  that  is  subjected  to  concurrent  dynamic 
pulses,  one  from  each  side,  for  a  maximally  dissipative  motion  is  examined  in  Lin  and  Pence 
(1991),  where  the  total  energy  loss  is  contrasted  to  that  which  would  occur  if  the  two  pulses  were 
not  concurrent  It  is  found  that  the  concurrent  pulse  encounter  suffers  the  greater  energy  loss  in  the 
event  that  both  incoming  pulses  ate  of  the  same  sign  (with  respect  to  strain),  whereas  the  concurrent 
pulse  encounter  suffers  the  lesser  energy  loss  in  the  event  that  the  incoming  pulses  are  of  opposite 
sign. 

The  ultimate  damping  capacity  for  these  dynamical  processes  is  dependent  upon  the  large-time  dis¬ 
sipative  behavior  of  the  ensuing  shear  wave  reverberations.  The  essential  obstacle  to  an  explicit 
analysis  of  this  problem  is  that  each  pulse-phase  boundary  encounter  spawns  both  a  reflected  pulse 
and  a  transmitted  pulse,  thus  giving  rise  to  distinct  pulse  generations,  involving  a  geometric  in¬ 
crease  in  the  number  of  pulses  which  must  be  treated  in  each  generation.  Faced  with  this  difficulty, 
we  developed  (Lin  and  Pence,  1992a)  an  analytical  framework  which  abandons  the  notion  of  track¬ 
ing  the  space-time  trajectory  of  each  individual  pulse,  and  instead  monitors  the  ongoing  change  in 
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the  stored  dynamical  field  energy  as  a  function  of  initiating  pulse  magnitude  and  reverberation  gen¬ 
eration  number.  Numerical  simulations  udliziog  this  scheme  have  enabled  us  to  examine  the  ques¬ 
tion  raised  above  regarding  the  possibility  of  completely  damping  out  an  initial  dynamical  pertur¬ 
bation.  We  establish  that  both  the  maximally  dissip^ve  motion,  and  the  maximal  energy  loss  mo¬ 
tion,  dissipate  exactly  the  amount  of  energy  needed  for  the  dynamical  system  to  setde  into  the  new 
energy  minimal  equilibrium  ground  state  associated  with  the  final  value  of  load.  The  numerical 
simulations  indicate  how  the  energy  loss  proceeds  by  reverberation  generation  number,  and  dem¬ 
onstrates  that  the  energy  remaining  at  each  generation  finely  distributes  itself  over  a  geometrically 
increasing  number  of  ever  smaller  pulses.  These  results  are  established  numerically  for  arbitrary 
sized  initiating  pulse;  we  were  also  able  to  obtain  the  same  results  analytically  for  the  special  case 
of  a  vanishingly  small  dynamical  perturbation.  This  lattra*  problem  is  rendered  meaning^  by  scal¬ 
ing  the  initiating  pulse  magnitude  with  respect  to  the  transformation  strain  before  extracting  the 
limit  In  fact  it  was  found  that  in  this  limit  each  generation  dissipates  exactly  one  half  of  the  ex¬ 
cess  energy  remaining  within  the  system.  This  intriguing  result  can  be  contrasted,  for  example, 
vdth  the  optimal  result  for  a  maximally  dissipative  motion  in  which  70.5%  of  the  energy  is  dissi¬ 
pated  within  the  first  generation.  These  results  are  apparently  the  first  of  any  kind  dealing  with  the 
large  time  dissipative  behavior  of  systems  governed  by  an  elastic  dynamics  supplemented  with  a 
local  criterion  for  phase  boundary  mobility.  As  such  they  form  a  useful  counterpart  to  the  results, 
given  recently  in  Ball  et  al.  (1991),  for  the  lai^e  time  dissipative  behavior  of  systems  governed  by 
field  equations  with  an  explicit  viscous  dissipation.  It  is  to  be  noted  that  the  dissipative  mechanism 
is  field  resident  in  Ball  et  al.  (1991)  and  hence  global  in  action,  whereas  in  our  case  the  dissipative 
mechanism  is  inherently  local. 

The  problems  discussed  in  Lin  and  Pence  (1991, 1992a)  have  focused  on  the  maximum  dissipation 
rate  criterion  for  determining  the  mobility  of  a  propagating  phase  boundary.  This  has  clarified  the 
ultimate  damping  capacity  of  these  materials  and,  in  particular,  indicates  how  dynamical  phase 
transition  processes  can  settle  down  into  equilibrium  states  governed  by  considerations  based  upon 
energy  minimization.  Phase  boundary  mobility  can  also  be  regarded  as  governed  by  a  deterministic 
kinetic  relation  between  the  phase  boundary  velocity  and  the  local  interface  driving  traction  (since 
the  latter  is  the  dissipation  rate  divided  by  the  phase  boundary  velocity,  this  is  also  equivalent  to  a 
relation  between  the  phase  boundary  velocity  and  the  dissipation  rate).  Noncontinuum  theories  (at¬ 
omistic  and  statistical  mechanical)  provide  guidance  in  the  formulation  of  such  relations.  The  prop¬ 
er  thermodynamical  framework  for  the  incorporation  of  a  kinetic  relation  has  been  clarified  by 
Gurtin  and  Struthers  (1990)  and  the  stability  of  kinetically  driven  phase  boundary  motion  has  been 
examined  by  Fried  (1991a,b).  Abeyaratne  and  Knowles  have  examined  the  role  of  a  kinetic  relation 
in  determining  phase  boundar7  motion  in  both  the  quasi-static  (1988a)  and  the  fully  dynamic 
(1991)  setting. 

To  the  extent  that  the  maximum  dissipation  rate  criteria  provides  a  relation  between  dissipation  and 
phase  boundary  velocity  (Pence  1991b)  it  can  be  regarded  as  a  kinetic  relation,  although  it  is  our 
belief  that  the  maximum  dissipation  rate  criterion  is  more  properly  viewed  as  providing  the  limiting 
behavior  for  dissipative  behavior  from  among  all  possible  kinetic  relations.  In  view  of  our  present 
understanding  of  this  limiting  behavior  as  given  in  Pence  (1991a,b)  and  Lin  and  Pence  (1992a), 
and  because  of  the  greater  generality  of  the  kinetic  relation  point  of  view,  we  have  begun  to  inves¬ 
tigate  the  extent  to  which  this  limiting  dissipative  behavior  is  attainable  from  kinetically  driven 
phase  boundary  propagation.  To  this  end,  in  Lin  and  Pence  (1992b)  we  have  reexamined  the  con¬ 
current  pulse  problem  for  the  case  in  which  phase  boundary  mobility  is  governed  by  a  kinetic  re¬ 
lation  instead  of  a  maximum  dissipation  condition.  Specifically,  we  have  considered  a  linear  rela- 
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tions  =  Kf  where  sand  fare  the  phase  boundary  velocity  and  the  phase  boundary  driving  traction, 
and  K  >  0  is  a  linear  coupling  constant  The  cases  k  -4  «>  and  k  -4  0  then  correspond  to  an  energy 
conserving  motion  and  a  rigidly  itmnobile  phase  boundary  respectively.  As  in  Lin  and  Pence 
(1991)  the  question  of  interest  is  to  determine  whether  the  greater  energy  loss  occurs  when  the  puls¬ 
es  act  concurrently  or  separately.  If  the  incoming  pulses  are  of  opposite  sign,  we  again  tind  that  the 
greater  energy  loss  occurs  when  the  pulses  act  separately.  However  if  the  incoming  pulses  are  of 
the  same  sign,  then  the  greater  energy  loss  is  experienced  by  the  concurrent  interaction  only  if  each 
of  the  incoming  pulses  is  sufficiently  small.  If  the  incoming  pulses  are  of  the  same  sign,  and  one 
or  the  other  is  sufficiently  large,  then  the  greater  energy  loss  takes  place  when  the  pulses  act  sepa¬ 
rately.  These  results  are  established  by  a  combination  analytical/numerical  analysis  which  allows 
us  to  calculate  the  threshold  values  on  pulse  size  which  lead  to  these  different  behaviors.  Work  is 
now  continuing  in  the  analysis  of  the  large-time  behavior  of  these  kinetically  driven  ringing  pro¬ 
cesses  and  it  is  clear  that  at  least  one  more  major  archival  paper  will  result  from  our  progress  to 
date. 

This  section  summarizes  results  that  are  reported  in  the  following  four  publications: 

Lin,  J.  and  T.  J.  Pence  (1991).  Energy  dissipation  in  an  elastic  material  containing  a  mobile  phase  boundary 
subjected  to  concurrent  dynamic  pulses.  Transactions  of  the  Ninth  Army  Conference  on  Applied  Matii- 
ematics  and  Computing.  F.  Dressel  ed.,  437-450. 

Lin,  J.  and  T.  J.  Pence  ( 1992a).  On  the  dissipation  due  to  wave  ringing  in  non-elliptic  elastic  materials,  (ac¬ 
cepted  for  publication  in  /.  Nonlinear  ScL). 

Lin,  J.  and  T.  J.  Pence  (1992b),  Kinetically  driven  elastic  phase  boundary  motion  activated  by  concunem 
dynamic  pulses,  to  appear  in  Transactions  of  the  Tenth  Army  Conference  on  Applied  Mathematics  and 
Computing.  F.  Dressel  ed. 

Pence,  T.J.  (1991b).  On  the  encounter  of  an  acoustic  shear  pulse  with  a  phase  boundary  in  an  elastic  mate¬ 
rial:  energy  and  dissipation.  J.  Elasticity  26, 95-146. 

In  addition,  another  manuscript  is  in  preparation  which  we  intend  to  submit  to  the  Journal  of  the 
Mechanics  and  Physics  of  Solids  in  late  Spring  1993. 

5.  Phase  Boundary  Creation  Due  to  Impact 

In  our  work  that  we  have  outlined  in  Section  4,  phase  transformations  were  due  to  the  motion  of  a 
pre-existing  phase  boundary.  Equally  important  is  the  generation  of  new  phase  boundaries.  To  clar¬ 
ify  this  aspect  of  phase  transformation  we  have  examined  impact  induced  phase  transitions  (Davi¬ 
son  and  Grahm  1979;  Asay  and  Kerley  1987).  Immediate  attention  was  restricted  to  models  for  iso¬ 
thermal  phase  transitions.  The  work  of  Slemrod  (1983, 1984)  and  Hattori  (1986a,b)  indicates  how 
these  isothermal  models  serve  to  clarify  the  overall  analytical  issues  before  treating  more  general¬ 
ized  models  (to  be  discussed  in  Section  6).  Even  in  this  preliminary  setting,  major  challenges  must 
be  confronted  due  to  the  possibilities  for  loss  of  hyperbolicity  in  the  governing  partial  differential 
equations  (Liu  1975;  Schaeffer  and  Shearer  1987;  Shearer,  Schaeffer,  Marchesin  and  Paes-Leme 
1987;  Issacson,  Marchesin,  Plohr  and  Temple  1988).  In  view  of  this  challenge,  it  is  common  to 
adopt  restrictions  on  the  size  of  the  variation  of  the  initial  data  (e.g.  Hattori  1986a).  However,  mo¬ 
tivated  by  situations  in  which  the  impact  velocity  may  be  large,  it  was  our  desire  to  obtain  a  meth¬ 
odology  that  does  not  hinge  upon  any  such  requirement  We  achieved  this  goal  for  the  fundamental 
Riemann  problem  by  focusing  on  separate  candidate  solutions  for  both  impactor  and  target  that  sat- 
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isfy  the  initial  conditions  of  the  problem.  These  solutions  are  then  matched  in  a  parameter  plane  of 
contact  interface  stress  a*  vs.  contact  interface  velocity  v*.  In  view  of  the  issues  involving  dynam¬ 
ical  nonuniqueness  mentioned  previously,  numerous  such  matches  may  be  possible,  all  of  which 
are  either  dissipative  or  else  energy  conserving.  In  a  space-time  diagram  (the  (x,t)-plane)  the  solu¬ 
tion  to  the  Riemaim  problem  partitions  itself  into  sectors,  and  in  each  sector  the  dynamical  fields 
may  either  involve  constant  strain  and  velocity  (a  C-sector),  or  else  a  continuous  compressive  wave 
(a  W-sector).  The  interface  between  sectors  may  be  continuous,  a  conventional  shock  (S)  or  else  a 
phase  boundary  (P).  We  first  considered  the  impact  between  a  target  material  that  admitted  phase 
transformations,  and  an  impactor  with  a  more  conventional  (albeit  nonlinear)  response  (Pence 
1990a,  1992a).  We  considered  a  range  of  impact  velocities  v  that  allow  a  variety  of  complicated 
wave  patterns  involving  phase  transitions  in  Ae  target  material.  In  this  regard  we  established  the 
existence  of  various  threshold  values  of  impact  velocity  V|  <  V2  <  V5  <  vj.)  (adopting  the  notation 
of  Pence  1992a)  such  that: 

(i)  For  0  <  V  <  V  j  there  is  a  single  solution  with  wave  pattern  (in  the  target)  of  type  C- W-C,  and,  in  view  of 
the  absence  of  shocks  and  phase  transitions,  this  solution  is  energy  conserving. 

(ii)  For  Vj  <  V  <  V2  there  are  exactly  two  solutions,  one  of  type  C-W-C  and  the  other  of  type 

C-W-C-P  (starting  from  the  x-axis  and  proceeding  to  the  t-axis).  Here  the  latter  differs  from  the  former 
(Mily  by  the  presence  of  a  stationary  (and  hence  non-dissipative)  phase  boundary  on  die  interface  of  contact 
(the  t-axis  in  a  Lagrangian  frame).  Both  such  solutions  are  energy  conserving. 

(iii)  For  V2  <  v  <  V5  there  is  a  one-parameter  family  of  solutions,  each  of  which  is  either  of  type 
C-W-C-P  or  C-W-P  or  C-W-P-C  or  C-W-C-P-C  or  C-P-C.  The  bounding  solutions  of  this  family  may  be 
energy  conserving,  while  the  interior  solutions  are  strictly  dissipative.  Thus,  at  most  2  of  these  solutions 
are  energy  conserving.  In  addition,  there  may  or  may  not  exist  a  single  (energy  conserving)  solution  of 
type  C-W-C. 

(iv)  For  V5  <  V  <  Vfj  there  is  exactly  one  solution  with  wave  pattern  of  type  C-P-C  and  this  solution  dissipates 
energy. 

The  impact  velocity  regime  v  >  v^  involves  con^licated  (single-phase)  wave  patterns  in  the  im¬ 
pactor,  but  no  additional  phase  transitions  in  the  target,  and  so  the  analysis  was  not  pursued  into 
this  regime.  The  far  more  significant  velocity  regime  (iii)  of  V2  <  v  <  V5  was,  however,  analyzed 
at  great  detail.  It  was  found  that  the  (o*,v’'')-plane  provides  a  natural  basis  for  investigating  the  ef¬ 
fect  of  different  solution  selection  criteria.  In  Pence  (1992a)  the  maximal  dissipation  selection  cri¬ 
teria  was  examined  by  determining  the  qualitative  features  of  the  level  curves  of  dissipation  in  the 
(<j*,v*)-plane.  It  was  shown  that  there  exists  another  threshold  velocity  V5  obeying  V2  <  V5  <  V5 
such  that: 

(iiia)  For  V2  <  v  <  V5  maximally  dissipative  solutions  are  restricted  to  the  interior  of  the  one-parameter  fam¬ 
ily  and  necessarily  have  wave  pattern  C-W-C-P-C.  Such  solutions  may  or  may  not  be  unique.  However, 
certain  special  forms  of  the  stress-response  frmction  for  the  target  material  wiU  ensure  uniqueness. 

(iiib)  For  V5  <  v  <  V5  there  is  exactly  one  maximally  dissipative  solution,  and  it  is  one  of  the  bounding  so¬ 
lutions  of  the  one-parameter  family.  The  resulting  wave  pattern  is  of  type  C-W-P-C. 

The  particular  special  forms  of  the  stress-response  function  which  were  shown  to  ensure  unique¬ 
ness  in  (iiia)  correspond  to  convexity  and  spacing  conditions  on  the  global  well  structure  of  the 
elastic  energy  density. 

The  extent  to  which  particular  dynamical  states  can  be  generated  (or  controlled)  can  be  formulated 
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as  an  inverse  Riemann  problem  (Pence  1990b),  Specifically,  one  seeks  to  determine  whether  a  giv¬ 
en  dynamical  wave  pattern  in  the  target  involving  a  moving  phase  boundary,  can  be  generated  by 
judicious  choice  of  impactor  material  and  impact  velocity.  Such  information  could  be  useful  for 
example  in  designing  explosive  cladding  processing/synthesis  operations  for  adaptive  composites 
involving  shape  memory  alloy  surface  structures.  The  particular  dynamical  states  that  are  so  attain¬ 
able  were  determined  in  terms  of  the  (a*,v*)-plane,  and  an  algorithm  was  given  for  the  determina¬ 
tion  of  the  associated  impactor  material  and  imp^t  velocity.  Although  the  particular  study  was 
limited  to  dynamical  phase  transformations  governed  by  the  maximum  dissipation  criterion,  the 
same  methodology  could  easily  be  applied  to  other  selection  criteria  such  as  those  involving  a  ki¬ 
netic  relation. 

Impact  between  two  materials,  each  of  which  is  capable  of  phase  transformations,  is  the  focus  of 
Pence  (1992b);  the  (a*,v*)-plane  remains  a  useful  organizing  principle.  In  this  case,  if  phase  trans¬ 
formations  take  place  simultaneously  in  the  impactor  and  target,  then  certain  ranges  of  impact  ve¬ 
locity  lead  to  a  two-parameter  family  of  dynamical  state  solutions  prior  to  the  invocation  of  a  dy¬ 
namical  selection  criterion.  Properties  of  the  level  curves  of  dissipation  were  obtained  for  this  prob¬ 
lem,  although  their  characterization  is  neither  as  simple  nor  as  complete  as  that  for  the  problem 
involving  only  a  single  phase  transforming  material  ^ence  1992a).  It  was,  however,  established 
that  the  maximally  dissipative  selection  criterion  may  lead  to  a  symmetry  breaking  bifurcation  phe¬ 
nomena  similar  to  that  which  occurs  in  Euler  column  buckling.  Recall  in  Euler  buckling  that  pairs 
of  “equally  likely”  asymmetric  buckled  states  (i.e.  right  buckled  &  left  buckled)  bifurcate  from  the 
symmetric  (unbuckled)  state  at  a  particular  value  of  load.  In  an  analogous  fashion,  the  symmetric 
impact  of  two  identical  phase  transforming  materials  may  give  symmetric  wave  patterns  for  low 
values  of  impact  velocity,  while  at  some  threshold  value  of  impact  velocity,  the  solutions  may  bi¬ 
furcate  into  a  pair  of  “equally  likely”  asymmetric  dynamical  wave  patterns  in  the  two  materials. 

This  section  summarizes  results  that  arc  reported  in  the  following  four  publications: 

Pence,  T.J.  (1990a).  Hiase  transitions  and  maximally  dissipative  dynamic  solutions  in  the  Riemaim  problem 
for  impact  Transactions  of  the  Eighth  Arniv  Confermce  on  Applied  Mathematics  and  Computing.  F. 
Dresseled.  817-832. 

Pence,  T.J.  (1990b).  An  inverse  Riemann  problem  for  impact  involving  phase  transitions.  Transactions  of 
the  Eighth  Army  Conference  on  Applied  Mathematics  and  Computing.  F.  Dressel  ed.  833-846. 

Pence,  T.J.  (1992a).  On  the  mechanical  dissipation  of  solutions  to  the  Riemann  problem  for  impact  involv¬ 
ing  a  two-phase  elastic  material.  Arch.  Ration.  Mech.  &  Anal.  117, 1-52. 

Pence,  T.J.  (1992b).  The  dissipation  topography  associated  with  solutions  to  a  Riemann  problem  involving 
elastic  materials  undergoing  phase  transitions,  to  r^pear  in  Shock  Induced  Transitions  and  Phase 
Structures  in  General  Media,  R.  Fosdick,  E.  Dunn  &  M.  Slemrod  eds.,  IMA  volume  series,  Springer- 
Verlag. 


6.  Thermal  Effects  Upon  Phase  Boundary  Mobility 

In  order  to  predict  the  temperature  change  associated  with  a  mechanically  activated  phase  transi¬ 
tion  we  have  extended  the  framework  of  Pence  (1991a,b)  in  order  to  encompass  the  role  of  thermal 
activation  or  inhibition.  The  utility  of  considering  such  extended  theories  is  demonstrated  by  the 
work  of  Slemrod  (1983, 1984)  and  Hattori  (1986a,b).  To  this  end,  we  have  developed  a  model  that 
is  thermodynamically  separable  in  that  it  can  be  embedded  into  a  fully  thermodynamically  consis¬ 
tent  theory  in  which  the  internal  energy  e  as  a  function  of  both  shear  strain  y  and  specific  entropy 
T]  additively  decouples  into  a  separate  functions  of  each  (this  idea  dates  back  to  Courant  and 
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Friedrichs  1948).  Mechanically  activated  motion  is  isothermal  within  such  materials.  In  order  to 
extend  our  modeling  capability  to  account  for  thermomechanical  coupling,  we  have  extended  the 
thermodynamical  framework  to  a  wider  class  of  materials  which  encompass  the  material  treated  in 
Pence  (1991a,b)  and  Lin  and  Pence  (1991, 1992a,b)  as  a  special  case.  The  research  described  in 
this  section  is  the  basis  of  graduate  student  Ralph  Worthington’s  Ph.D.  research. 

We  consider  a  material  in  which  each  phase,  say  indexed  by  i,  has  internal  energy  e,  stress  x  and 
entropy  7)  as  functions  of  shear  strain  y  and  temperature  6  given  by 

Ei  =  ^lij[Y-ai]^  +  bi0  +  bi ,  =  iijy-ail  -kiS,  tIj  =  bjlnO  +  kiY+ki , 

where  a^,  b^,  bj,  kj  and  kiare  constitutive  parameters  which  determine  the  spacing  and  curva¬ 
ture  of  the  individual  phase  energy  wells,  as  well  as  the  thermomechanical  coupling.  In  particular, 
kj  is  a  thermomechanical  coupling  parameter  since  any  phase  for  which  k^  =  0  is  governed  by  the 
purely  mechanical  theory  treated  in  Pence  (1991a,b)  and  Lin  and  Pence  (1991, 1992a,b).  The  equa¬ 
tions  of  compatibility,  force  balance  and  energy  balance  give  a  system  of  partial  differential  equa¬ 
tions  similar  to  that  obtained  by  Niezgodka  and  Sprekels  (1988)  in  their  studies  of  phase  transfor¬ 
mations,  albeit  of  different  form  since  the  energy  potential  is  not  of  Landau-Devonshire  type.  More 
signiflcandy  it  does  not  allow  field  resident  viscous  dissipation,  and  there  is  also  the  assumption 
that  heat  conduction  is  negligible  compared  to  the  other  process  variables,  thus  having  the  advan¬ 
tage  of  yielding  a  model  with  distinct  phase  interfaces.  This  model  is  in  fact  qualitatively  similar 
to  the  general  class  of  models  described  by  James  (1983),  in  addition  it  also  incorporates  the  cou¬ 
pling  parameter  kj  which  enables  us  to  investigate  how  the  thermomechanical  theory  unfolds  from 
a  purely  mechanical  theory  whose  behavior  is  well  known. 

We  have  examined  the  encounter  of  a  pulse  with  a  phase  boundary  in  this  thermomechanical 
framework  for  the  case  in  which  k]  =  0,  and  k2  >  0  corresponding  to  a  thermally  active  second 

phase.  Consider  for  the  moment  a  thermally  inactive  second  phase  (both  kj  =  0,  and  £3  =  0). 
Then  the  pulse-phase  boundary  encounter  problem  requires  the  consideration  of  two  regions  in  a 
space-time  plane:  a  region  corresponding  to  the  transmitted  pulse  and  a  region  corresponding  to  an 
interaction  zone  (bounded  by  the  incoming  pulse,  the  phase  boundary,  and  the  reflected  pulse).  By 
means  of  characteristic  theory,  four  fundamental  (linear)  relations  between  the  strain  and  velocity 
in  these  two  regions  can  then  be  obtained  (Pence  1991a).  When  these  equations  are  solved  in  terms 
of  the  strain  magnitude  of  the  incoming  pulse  Ay,  the  phase  boundary  speed  is  not  determined  and 
hence  is  obtained  by  an  external  mobility  criterion  of  the  type  discussed  in  Truskinovsky  (1990), 
Pence  (1991a,b),  and  Abeyaratne  and  Knowles  (1991).  Let  us  now  return  to  the  case  of  a  thermally 
active  second  phase  k2  >  0  and  assume  that  the  pulse  is  incident  through  the  thermally  inactive  fu^t 
phase  and  that  the  phase  boundary  is  driven  into  the  thermally  active  second  phase.  One  then  finds 
that  the  interaction  zone  now  splits  in  two  along  a  (new)  characteristic  associated  with  the  temper¬ 
ature  variable.  There  are  then  nine  fundamental  field  quantities  (strain,  velocity  and  temperature  in 
each  of  the  three  regions).  By  means  of  characteristic  theory,  we  have  obtained  nine  equations  re¬ 
lating  these  field  quantities  to  the  disturbance  magnitude  Ay  and  the  phase  boundary  speed  s .  These 
equations  are  not  linear  in  the  temperatures.  We  have  been  able  to  reduce  these  equations  to  a  single 
highly  nonlinear  algebraic  equation  for  the  tentperatureOj  in  the  transmitted  wave  adjacent  to  the 


phase  boundaiy: 
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r^(0)  -^(0)  1 

wheref(0)  =  2<D(0)+«&(O)ln[^i^^-j^^Jand«I>(0)  = 

material  mass  density  and  0  is  the  ambient  temperature  in  the  material  before  the  encounter.  The 
results  given  in  Pence  (1991a)  with  Cj  =  Jii/p  must  be  retrievable  in  the  limit  k2  ->  0.  Prelimi¬ 
nary  indications  are  that  this  is  indeed  the  case,  for  example  we  find  for  small  £2  that  the  above 
equation  yields 


+  0 


.  Here  pis  the 
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where  r  =  (£7lhP“lii)2AY+sJp^(Yi-Y2)  and  ^  =  p^J^+J^^^J^-sJwethen 
find  that  a  long  calculation  gives,  for  example,  that  the  velocity  in  the  transmitted  wave  is 
Vj  =  Vj  +  Vj^  where  v™ is  the  velocity  for  the  problem  in  which  k2  =  0: 


,,mech  _ 
Vt 


and  Vj^,  which  is  order  k^,  is  given  by 


ri^rS'ifi ,  ,V  1 

U,4  J[5JpI  2jir,  )  2j^,_ 


+  o(^). 


As  expected,  the  leading  order  term  in  the  above  result  coincides  with  that  given  in  Pence  (1991a) 
once  allowance  is  made  for  the  pulse  being  incident  through  phase- 1.  We  anticipate  that  all  of  the 
nine  field  quantities  will  similarly  yield  to  a  formal  asymptotic  series  where  the  leading  order  term 
retrieves  the  results  presented  in  Pence  (1991a).  These  results  will  enable  us  to  precisely  gauge  the 
extent  to  which  purely  mechanical  theories  of  phase  boundaiy  motion  can  be  unfolded  into  fully 
thermomechanical  theories.  An  assessment  of  the  qualitative  robusmess  of  this  unfolding  process 
then  demands  an  inquiry  as  to  whether  or  not  completely  new  solutions  appear  during  this  unfold- 
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ing  process.  These  questions,  and  related  issues,  are  currently  under  active  investigation. 

The  results  presented  in  this  section  summarize  the  results  of  collaboration  with  Ph.D.  graduate  stu¬ 
dent  Ralph  Worthington.  Ultimately,  in  addidon  to  Mr.  Worthington’s  Ph.D.  thesis,  we  estimate 
that  at  least  two  archival  papers  will  result,  which  at  present  we  are  targeting  for  the  SIAM  J.  of 
Applied  Math,  and  for  Continuum  Mechanics  and  Thermodynamics. 

7.  Thermal  and  Mechanical  Evolution  of  Phase  Variants 

All  of  the  results  discussed  thus  far  are  devoted  to  the  modeling  of  phase  boundary  motion  at  the 
micromechanical  level.  The  leap  to  a  macromodel  requires  the  consolidation  of  the  various  effects, 
describable  individually  at  the  microlevel,  into  a  constitutive  theory  which  accounts  for  their  cu¬ 
mulative  effect  at  a  level  of  detail  appropriate  for  the  model’s  application  area.  To  this  end,  we  have 
been  collaborating  with  an  industrial  counterpart  (Y.  Ivshin  of  Johnson  Controls)  in  order  to  bridge 
this  gap  in  a  fashion  that  incorporates  the  ever  increasing  knowledge  base  obtained  from  the  type 
of  detailed  micromechanical  modeling  discussed  thus  far. 

Both  the  shape  memory  effect  and  pseudo-elastic  effect  are  examples  of  hysteretic  response  in  the 
austenite  A  to  martensite  M  transition.  For  a  given  sh£q)e  memory  alloy  active  element,  this  re¬ 
sponse  can  be  characterized  experimentally  for  complete  transition  paths  (pure  A  to  pure  M  and 
vice-versa)  by  various  standard  techniques  (calorimetry,  x-ray  difb^tion,  resistivity,  load-cell). 
The  pure  A  to  pure  M  transition  will  thus  give  rise  to  a  response  surface  in  the  relevant  parameter 
space  of  design  variables  (e.g.  temperature,  stress,  displacement),  and  the  pure  M  to  pure  A  tran¬ 
sition  will  give  rise  to  a  similar,  but  distinct,  surface.  Together  these  two  surfaces  form  an  envelope 
which  bounds  the  system  response  for  all  possible  transition  paths.  Its  significance  derives  from 
the  fact  that,  within  a  precise  sensing  and  control  application,  an  active  element  will  invariably  un¬ 
dergo  partial  transitions  before  arresting  itself  and  reversing  its  response,  thus  generating  a  com¬ 
plicated  trajectory  within  the  bounding  envelopes.  Algorithms  for  predicting  these  trajectories 
from  the  knowledge  of  the  complete  transition  behavior  are  of  fundamental  concern,  as  it  is  im¬ 
practical  to  experimentally  compile  a  complete  data  base  description  of  all  arrested/reversed  tran¬ 
sition  behavior  for  a  given  shape  memory  alloy  active  element  (i.e.  to  experimentally  determine  the  tra¬ 
jectories  of  all  space  filling  curves  within  the  bounding  envelopes).  Analytical  characterization  of  these  tra¬ 
jectories  is  also,  of  course,  central  to  a  continuum  based  constitutive  description  of  such  material 
behavior. 

Various  algorithms  have  been  proposed  for  the  modeling  of  shape  memory  material  constitutive 
response  (e.g.  Muller  1979;  Muller  and  Wilmanski  1980;  Liang  and  Rogers  1990;  Zhang,  Rogers 
and  Liang  1991),  but  in  general  the  precise  characterization  of  arrested/reversed  transition  behavior 
in  terms  of  complete  transition  behavior  has  not  been  the  central  concern.  When  it  has  been,  issues 
related  to  proper  qualitative  algorithmic  behavior  have  not  been  completely  addressed,  thus  possi¬ 
bly  leading  to  models  which,  for  example,  allow  certain  response  trajectories  involving  arrested/ 
reversed  transitions  to  exit  the  bounding  envelopes.  In  Ivshin  and  Pence  (1992a)  we  have  given  ex¬ 
amples  of  how  seemingly  reasonable  hysteresis  algorithms  can  give  rise  to  such  cases  of  envelope 
piercing.  To  address  such  issues  we  have  developed  a  model  for  rate  independent  hysteresis,  gov¬ 
erned  by  differential  equations  of  Duhem  Madelung  type  (Visintin  1988),  which  is  particularly  ap¬ 
plicable  to  tenq)erature  induced  phase  transitions  which  occur  in  shape  memory  alloys.  Other  phys¬ 
ical  systems  exhibiting  hysteresis  have  been  modeled  within  a  similar  framework,  notably  ferro- 
magnetically  soft  materials  (Coleman  and  Hodgdon  1986, 1987).  In  its  simplest  form,  as  examined 


14 


in  Ivshin  and  Pence  (1992a),  the  model  treats  the  evolution  of  a  relative  phase  fraction  variable  (e.g. 
austenite  fraction)  due  to  a  prescribed  temperature  history.  Micromechanical  considerations,  nota¬ 
bly  that  partial  transitions  proceed  at  a  pace  that  is  proportional  to  the  phase  fraction  of  the  parent 
phase  (the  phase  A  or  M  that  is  being  depleted),  are  ensured  by  the  form  of  the  governing  relations.  This 
leads  to  a  correspondence  between  experimentally  determinable  transition  envelopes  for  complete 
transitions,  and  thermodynamic  phase  fraction  transition  rate  functions. 

We  believe  that  this  model  extends  previous  modeling  efrbrts  in  three  key  respects.  First,  connec¬ 
tion  to  experiment  follows  directly  from  the  dependence  upon  a  hysteresis  envelope  for  complete 
transitions,  and  envelope  restrictions  have  been  obtained  which  ensure  proper  monotonicity,  con¬ 
tainment  and  orientation  requirements  upon  partial  transition  paths.  Second,  we  have  identified  un¬ 
derlying  mathematical  structures  which  rigorously  characterize  the  hysteresis  response.  These 
structures  include:  attracting  limit  cycles  for  finite  cycling,  and  attracting  limit  points  for  infinites¬ 
imal  cycling  (drift).  These  structures  are  anticipated  to  allow  for  the  rational  analysis  of  seemingly 
chaotic  response  trajectories.  Third,  although  many  models  for  the  shape  memory  affect  are  capa¬ 
ble  of  capturing  isothermal  stress-strain  hysteresis,  the  generalization  to  adiabatic  and  convective 
heating  has  not  received  appropriate  attention,  especially  since  the  latter  is  critical  for  device  sim¬ 
ulation.  We  have  successfully  simulated  isothermal,  adiabatic,  and  convective  stress-strain  behav¬ 
ior  in  a  simple  prototype  model.  In  particular,  tiie  general  convective  stress-strain  behavior  suc¬ 
cessfully  approaches  adiabatic  behavior  for  low  convectivity/fast  loading,  and  successfully  ^ 
proaches  isothermal  behavior  for  high  convectivity/slow  loading.  These  results  also  allow  for 
simulation  of  multiple  heating/cooling  mechanisms  acting  simultaneously  (e.g.  resistive  heating/ 
convective  cooling).  Since  the  prepublication  manuscript  which  details  this  successful  simulation 
of  convectivity  behavior  is  still  in  preparation,  we  present  a  short  summary  diagram  in  Figure  1. 

In  Ivshin  and  Pence  (1992b),  we  have  implemented  the  minimum  amount  of  extensions  to  Ivshin 
and  Pence  (1992a)  that  allow  for  a  micromechanically  consistent  modeling  of  the  two-way  mem¬ 
ory  effect  This  requires  the  consideration  of  two  martensite  variants:  M+  and  M-,  with  different 
transformation  strains.  The  underlying  phase  transitions  between  the  three  phase/variant  species  A, 
M-f  and  M-  are  then  triggered  both  mechanically  and  thermally  on  the  basis  of  changes  in  material 
stability  as  governed  by  an  enlarged  set  Duhem-Madelung  differential  equations.  In  Ivshin  and 
Pence  (1992b),  we  considered  a  simple  symmetric  three  layer  composite  system,  each  layer  of 
which  is  composed  of  the  same  shape  memory  material.  Unequal  initial  lengths  prior  to  bonding 
gives  rise  to  a  piecewise  constant  residual  stress  distribution  as  one  proceeds  through  the  layers. 
For  this  simple  system,  the  two  way  memory  effect  manifests  itself  in  terms  of  length  changes 
which  occur  upon  raising  and  lowering  the  system  temperature.  This  behavior  stems  from  the  abil¬ 
ity  of  the  model  to  account  for  asymmetry  in  transitions  between  the  three  phase/variant  species 
owing  to  the  nonuniform  nature  of  the  residual  stress  distribution.  A  complete  description  of  this 
system  required  the  consideration  of  eleven  variables:  the  conunon  system  temperature;  the  stress, 
strain  and  phase  fractions  of  A,  M-i-  and  M-  in  the  (synunetric)  outer  layers;  and  the  stress,  strain 
and  phase  fractions  of  A,  M-f  and  M-  in  the  middle  layer.  Ten  differential  equations  relating  these 
variables  follow  naturally  from  the  model,  which,  when  augmented  with  a  prescribed  applied  tem¬ 
perature  T(t),  gives  rise  to  a  well  defined  mathematical  system.  Numerical  simulations  fora  peri¬ 
odic  T(t)  exhibited  cyclic  length  changes  to  the  overall  system,  which  rapidly  collapsed  onto  a  limit 
cycle  when  graphed  against  temperature.  The  other  system  variables  also  exhibited  a  rapid  collapse 
onto  their  associated  limit  cycles. 
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This  section  summarizes  results  that  are  reported  in  the  following  two  publications: 

Ivshin,  Y.  and  T.  J.  Pence  (1992a).  A  constitutive  model  for  hysteietic  phase  transition  behavior,  (accepted 
for  publication  in  Ira.  J.  Eng.  Sci.). 

Ivshin,  Y.  and  T.  J.  Pence(1992b).  A  simple  mathematical  model  of  two-way  memory  effect,  (to  appear  in 
Proceedings  of  die  1992  Int.  Corf,  on  Martensitic  Tranrformations  (ICOMAT-W)). 

In  addition,  another  manuscript  is  substantially  complete  which  we  intend  to  submit  in  early  Spring 
1993  to  either  the  Int.  J.  Eng.  Sci,  or  to  J.  Intelligent  Materials  and  Structures. 

8.  Laminated  Structures  and  Thin  Films 

Within  the  setting  of  finite  deformation  elasticiQr,  the  ability  of  a  particular  material  model  to  admit 
the  abrupt  nucleation  of  phase  boundaries  and  other  microstructural  features  responsible  for  qual¬ 
itatively  different  macroscopic  response  is  dependent  upon  the  constitutive  energy  function  of  the 
material.  In  this  regard,  material  instabilities  at  the  macrolevel  originate  as  structural  instabilities 
at  the  microlevel.  The  associated  bifurcation  phenomena,  whose  behavior  may  be  well  understood 
in  monolithic  systems,  can  drastically  change  its  character  in  composite  systems.  In  highly  deform¬ 
able  thermoelastically  “passive’'  materials  (they  do  no  not  admit  phase  changes),  our  previous  work 
has  shown  how  the  buckling  of  layered  structures  leads  to  configuration  modification,  and  buckling 
load  reordering,  for  those  buckled  states  that  ate  natural  extensions  of  previously  known  solutions 
corresponding  to  an  isolated  layer.  For  the  isolated  layer  it  had  been  shown  (Sawyers  and  Rivlin 
1974)  that  there  are  two  families  of  buckled  solutions:  a  conventional  andsyirunetric  flexural  buck- 
Ung  mode  (reminiscent  of  classical  Euler  buckling),  and  a  symmetric  barreling  buckling  mode.  The  buck¬ 
ling  loads  associated  with  these  different  buckled  solutions  are  ordered: 

Tf  <  T?  <  ...<Tf  <  <...  .^  T.  ,<  T»  <  <  Tf, 

where  T^  and  T^  are  the  buckling  loads  associated  with  the  m^^  flexure  and  barreling  mode  re¬ 
spectively,  and  T^  is  associated  with  a  wrinkling  instability.  In  contrast,  for  a  three-ply  sandwich 
construction  we  find  certain  instances  in  which  T^  is  the  minimum  buckling  load  so  that  wrin¬ 
kling,  and  not  mode-1  flexure,  is  the  critical  instability  (Pence  and  Song  1991). 

These  buckling  results  have  immediate  consequences  for  the  optimal  design  of  laminated  struc¬ 
tures  with  respect  to  buckling  resistance.  This  was  demonstrated  in  Song  and  Pence  (1992)  for  the 
following  simple  pilot  problem:  Construct  a  symmetric  three-ply  layered  composite  of  given  di¬ 
mensions  from  two  neo-Hookean  materials,  A  and  B,  whose  combined  amount  is  equal  to  the  total 
volume.  Thus  either  AorB  may  be  used  for  the  central  ply.  The  optimal  construction  is  that  which 
gives  the  largest  critical  thrust.  We  found  that  the  optimal  design  was  dependent  on  the  aspect  ratio 
of  the  construction.  If  the  structure  is  short  in  the  direction  of  thrust,  then  the  optimal  construction 
was  one  in  which  the  stiffer  material  was  used  in  the  central  ply.  However,  if  the  structure  is  long 
in  the  direction  of  thrust,  then  the  other  construction  is  optimal.  Thus  there  is  a  transition  aspect 
ratio  in  this  simple  optimal  design  problem.  Numerical  routines  were  developed  to  determine  this 
aspect  ratio.  In  addition,  by  introducing  composite  construction  characterization  parameters  p  and 
a,  where  |3  is  the  ratio  of  the  ply  stiffnesses  and  a  is  the  volume  fraction  of  the  central  ply,  we  suc¬ 
cessfully  implemented  a  (P,a)-parameter  plane  characterization  of  the  optimal  design  criteria. 

Recently,  we  have  established  the  existence  of  completely  new  buckled  states  that  have  no  isolated 


16 


layer  counterpart  For  example,  it  appears  for  an  N-layered  structure  that  there  exists  N- 1  new  fam¬ 
ilies  of  buckled  solutions,  in  addition  to  the  2  previously  known  families  corresponding  to  flexure 
and  barreling  (we  have  verified  this  for  N=2,3  while  for  N>3  this  is,  at  present  our  working  hy¬ 
pothesis).  We  are  developing  appropriate  numerical  algorithms  for  unraveling  the  complicated  bi- 
fiircation  behavior  and  clarifying  the  properties  of  these  new  solutions.  In  particular,  asymptotic 
analysis  has  enabled  us  to  determine  precise  load  levels  which  lead  to  wrinkling  phenomena.  In 
addition  these  new  deformations  will  be  neither  pure  flexure  nor  pure  barreling,  and  we  are  devel¬ 
oping  methods  for  characterizing  this  mode-mixity  (the  importance  stems  from  the  fact  that  the  dif¬ 
ferent  mode  types  may  be  prone  to  different  delaminadon  risk).  Although  the  analysis  is  at  present 
limited  to  passive  nonlinear  elastic  materials,  we  are  here  building  toward  the  goal  of  analyzing 
multi-layers  that  involve  both  passive  and  thermoelasdcally  acdve  layers,  since  the  presence  of  the 
latter  offers  the  possibility  for  acdve  mechanical  tuning  of  an  overall  smart  stmcture.  The  results 
given  in  Ivshin  and  Pence  (1992b)  indicate  that  this  can  be  accomplished  in  thermoelasdcally  ac¬ 
dve  muld-ply  structures  for  simple  uniaxial  deformadons,  whereas  here  we  are  concerned  with 
more  complicated  buckling  phenomena  as  well  as  more  general  out  of  plane  displacement  modes. 

The  results  of  our  study  of  laminated  structures  are  applicable  to  thin  films  in  an  appropriate  lim- 
idng  case.  As  a  pracdcal  matter,  thermoelasdc  phase  transformadon  behavior  is  understood  far  bet¬ 
ter  for  materials  in  bulk  form  than  it  is  for  thin  film  structures.  Here  the  influences  low  dimensional 
constraint;  specialized  fabricadon  induced  variant  structures  (and  their  associated  residual  stress 
fields);  and  the  compeddon  between  translaminar  and  intralaminar  heat  flow  all  serve  to  compli¬ 
cate  the  structural  response.  We  have  recendy  begun  to  invesdgate  the  possibilides  for  modeling 
these  phenomena  (an  avenue  of  study  which  was  completely  unforeseen  in  1989  when  the  present 
study  was  proposed).  We  hope  to  pursue  this  invesdgadon  via  coUaboradon  with  materials  science 
colleague  David  S.  Grummon.  Recendy,  Prof.  Grummon  has  successfully  fabricated  NiTl  amor¬ 
phous  thin  films,  and  subsequendy  has  developed  annealing  methods  to  achieve  crystallogrq>hic 
structure  which  enable  thermoelasdc  austenit^martensite  transformadons  that  underlie  shape- 
memory  and  superelasdcity  effects.  This  is  the  radonale  for  our  joint  proposal  to  ARO,  which  was 
mendoned  earlier  in  this  report. 

This  section  summarizes  results  that  are  reported  in  the  following  publication: 

Song,  J.  and  T.  J.  Pence  ( 1992),  On  the  design  of  three-ply  nonlinearly  elastic  composite  plates  with  optimal 
buckling  performance.  Structural  Optimization  5  (1992)  45-54. 

In  addition,  another  manuscript  is  substantially  complete  which  we  intend  to  submit  in  late  winter 
1993  to  Int.  J.  Solids  and  Structures. 
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